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Abstract 

For the nonlinear Schrodinger equation (NLSE), in presence of disorder, exponentially localized 
stationary states are found. In the present work it is demonstrated analytically that the localization 
length is typically independent of the strength of the nonlinearity and is identical to the one found 
for the corresponding linear equation. The analysis makes use of the correspondence between the 
stationary NLSE and the Langevin equation as well as of the resulting Fokker-Planck equation. 
The calculations are performed for the "white noise" random potential and an exact expression 
for the exponential growth of the eigenstates is obtained analytically. It is argued that the main 
conclusions are robust. 

PACS numbers: 72.15.Rn, 42.25.Dd, 42.65.k 
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In this work we consider a simple problem of the one- dimensional Anderson localization 
l|, |2] for the nonlinear Schrodinger equation (NLSE) . This problem is relevant to 

experiments in nonlinear optics, for example disordered photonic lattices 6|], where Anderson 
localization was found in presence of nonlinear effects^aswell as experiments on Bose-Einstein 
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Condensates (BEC) in disordered optical lattices 
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disorder and nonlinear effects leads to new interesting physics 
^articular, the problem of spreading of wave packets and transmission are not simply related 
la, 113], in contrast with the linear case. In spite of the extensive research, many fundamental 



problems are still open, and, in particular, it is not clear whether in one dimension (ID) 
Anderson localization can survive the effects of nonlinearities. 

Herein we consider ID localization of stationary solutions of the NLSE in a random 
potential. The problem is described by the equation 



(1) 



where V{x) is a random 5-correlated potential with a Gaussian distribution, of zero mean 
and variance a^, such that 



{V{x)V{x)) = 2a^6{x - x) 



(2) 



where (. . .) denotes the average over realizations of the random potential. The variables are 
chosen in dimensionless units and the Planck constant is h = 1. For the linear case {(3 = 0) 
this model was studied extensively in the past jisl. [lol. The problem in question is Anderson 
localization of stationary solutions of Eq. with energies u 



ipixjt) = exp{—iujt)(f){x) 



where (f){x) is real. Substituting Eq. ([3]) in Eq. ([T]) one obtains the stationary NLSE 



ujcf){x) = -dl<j){x) + (3<f{x) + V{x)(P{x) 



(3) 



(4) 



It was established rigorously 3|, |20|, l21| that this equation has exponentially localized solu- 
tions of this type for a wide range conditions. It is instructive to notice that also in absence 
of disorder {V = 0) the nonlinear equation exhibits stationary localized states, and these are 



not simply related to the localized states in presence of disorder [13|, [IJ, l22j . The purpose of 



the present work is to find the localization length of these states. It will turn out that the 



localization length is typically not affected by the nonlinearity and is identical to the one 



of the linear problem [jS = 0). The approach d la Borland 23|, |2J] will be the basis of our 
analysis. This approach is reviewed clearly in detail in 25| and made rigorous in 261 ] . 

We will specifically calculate of solutions of Eq. (jlj) that are found for a certain 

u, with given boundary conditions at some point, for example 0(x = 0) and = 0), 
where prime means the derivative with respect to x. It will be shown that this quantity 
grows exponentially with the rate 

27= lim ^ >0, e = -, 5 

that is independent of (3, where ^ is the localization length. Note, it is different from 
the usually studied self averaging quantity 7^ = |^(ln0^(x)) = | limj-^oo ^'^'^J^^ . We will 
find that 7 is a smooth function of energy. Since the distribution of random potentials is 
translationally invariant, it is independent of the choice of the initial point as a; = 0. Like in 
the linear case, starting from a specific initial condition, will typically grow. For specific 
values of u at some point this function will start to decay, so that a normalized eigenfunction 



is found 



23l . |2J, |25|, |26| . The envelope of the wave function will grow exponentially if we start 



either from the right or from the left. The value of u results from the matching condition, 
so that an eigenfunction has some maximum and decays in both directions as required by 
the normalization condition. The exponential decay is an asymptotic property, while the 
matching is determined by the potential in the vicinity of the maximum. This observation 



25 



26l | is crucial for the validity of this approach and enables us to determine the exponential 



decay rate of states from the solution of the initial value problem (jlj). For the linear case 
these values of u form the point spectrum of the problem that is the entire spectrum of the 
linear problem. This approach can be followed also for the nonlinear problem, but contrary 
to the linear case, these stationary states do not provide a complete picture of the dynamics. 
Let us fix the stationary states in the vicinity of u will on the average (over realizations) 
decay with the localization length ^ defined by Eq. ([5]). As mentioned it will be found to be 
independent of /?, leading to the conclusion that this is a typical property of the localized 
eigenstates of Eq. (jlj). 

The calculation of {(f)'^{x)) will be performed by the analogy with the classical Langevin 
equation jisl. fioll. Therefore, here we are considering the x-coordinate as the formal time on 
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the half axis x = t E [0, oo), and Eq. (jlj) reduces to the Langevin equation 

4) + uj(j)- (3(f)^ -V{T)(j) = (6) 

with the 6 correlated Gaussian noise V{t). Now we introduce new variables u = cj) and 
V = <f) = ^ and a distribution function of these new variables is P = P{u,v,t). The 
dynamical process in the presence of the Gaussian 5-correlated noise is described by the 



distribution function that satisfies the Fokker-Planck equation: (FPE) (see, e.g., 271. l28l|) 



drP - [lou - (3u^]d^P + vduP - a'u'dtP = . (7) 



It is obtained from the Langevin equation ([6]), as can be seen in 



with [19], Eq. (6.14) there |29|. 
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28| and is consistent 



We are interested in the average quantum probability density (0^(x)) = {u^{t)), where 

{u^{t)) = j u^P{u,v,T)dudv . 

It is useful to obtain from the FPE a system of equations for the moments 

Mk,i = {u%') , (8) 

where k,l = 0,1,2,.... Substituting u'^v^ in the FPE and integrating over u and v, one 
obtains the following relation for ^ 

Mk,i = -lujMk+i,i-i + A;Mfc„i,,+i + /(/ - l)a^Mk+2,i-2 + l3lMk+3,i-i , (9) 

where Mk^i with negative indexes are assumed to vanish. We note that only terms with the 
same parity of A; + Z are coupled. Since we are interested in M2,o = (m^), we study only the 
case when this parity is even, namely k + 1 = 2n with n = 1,2, . . .. The sum of the indexes 
of the moments is 2n, except the last term (3lMj.j^3^i_i, where the sum is 2{n + 1). This leads 
to the infinite system of linear equations that can be written in the form 

M = >VM, (10) 

where M = (M2,o? ^0,2, ^4,0? ^3,i; • • •) = W is the corresponding matrix. The 
matrix elements yVk,i are determined by Eq. ([9]). The solutions of the system of linear 
equations (ITOl) are linear combinations of the eigenfunctions 

M^{t) = exp(Ar)M;,(0) , (11) 



where Ma(0) is the eigenvector of W corresponding to A. The growth rate of each moment, 
in particular M2,o = (m^), is determined by the eigenvalue with the largest real part Re A. 

For (3 = 0, Eq. (or (ITU]) ) has a closed form for each n. Therefore, the infinite matrix 
W is block diagonal and consists of the independent blocks An = An[{2n + 1) x [2n + 1)], 
and the characteristic polynomial reduces to a product of their determinants 

oo 

nciet(A„-A/„)=0, (12) 

n=l 

where J„ is an (2n+ 1) x [2n + 1) unit matrix. The problem of localization in the framework 
of Eq. ([9]) reduces to n = 1. The relevant characteristic polynomial det {Ai — A/i) reduces 
to a cubic equation 

A^ + 4cjA - 4(7^ = . (13) 

Cardano's method yields jsol 

X, = R+ + R_, A2,3 = ^ ± , (14) 



where R± 



2a' ± ./4a4 + 



. The growth rate is determined by the eigenvalue with 



27 

largest real part that will be denoted by A^. We conclude that asymptotically for large x 
the averaged wave function indeed grows exponentially as (0^(x)) ~ e'^'^^ and 27 = Am- Also 
the behavior of the higher blocks can be calculated. From the n-th block the behavior of 
the moment (0^"(x)) can be found. In the high energy limit one finds from Eq. iHM 

Am = Ai ^ — , uj ^ +00 , (15) 

UJ 



Am = Re A3 ~ 2-i/ |c<j| , 00^—00. (16) 



These limits can be found directly from Eq. 



T3l) . The solutions should be compared with 



the high energy asymptotics obtained in Ref 19|] (Eq. (10.12), p. 143), where 7^ = a'^/Auj, 
and 7s = \f\uj\ in limits u +00 and u —00, respectively. Note that Eq. (IT^ gives a 
simple expression for 7 for all values of the parameters, while the expression for 7^ is known 
only in the large u limit. Since 7 and 7^ result of different averages, these are not expected 
to be identical. 

Now, let us consider localization for /3 7^ 0. The eigenvalue equation 

W{(3)Mx{(3) = mMx{(3) , (17) 
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is obtained from Eq. ([9]), W is not block diagonal anymore. The (3 dependence results from 
the last term in Eq. ([9]). The (3 dependent terms couple the n-th block and (n + 1) block 
and are located above the {n + 1) block and to the right of the n-th block. Consequently, 
the (3 dependent terms do not affect the characteristic polynomial, as can be shown by 
elementary operations on determinants. It reduces to the one found for /? = 0, namely Eq. 
( IT2|) . Therefore, growth rates of all moments of do not depend on (3 and their values 
are equal to the ones of the linear problem for (3 = 0. This is correct in particular for Am, 
consequently 7 of Eq. ([5]) is identical to the value found in the linear case (/5 = 0). 

We demonstrated that in the presence of a random potential the stationary states of the 
NLSE are exponentially localized with the localization length that is found in the absence 
of the nonlinearity. This is in agreement with a heuristic argument that the effect of non- 
linearity is negligible where the wave function is small. We believe that the approach d 
la Borland can be extended to a rigorous treatment of the stationary states of the NLSE. 
Since this equation is nonlinear, the stationary states do not provide the complete or even 
essential description of the dynamics, starting from a given initial condition. The status of 
stability of these states with respect of small perturbations is not clear. The relation to the 
transmission problem is not obvious. Our results are consistent with the limit of vanishing 
flux in the transmission problem 31|. 

Another question that should be discussed is of the generality of the results. Assume 
that the power of ip in Eq. ([T]) differs from 3. Only the last term in Eq. ([9]) will be 
affected resulting in a different coupling between the blocks of the matrix W. But since 
these couplings are above the diagonal of the block diagonal matrix, they will not affect 
our conclusion that the characteristic polynomial which determines A is not affected by 
the nonlinearity. If the potential V{x) deviates from a white noise one, Kramers-Moyal 



coefficients 



27( 1 that are higher than the second one appear. But if no convergence problems 



of the Kramers-Moyal expansion are encountered, an equation like Eq. ([7]) with higher 
powers of u and f , combined with higher order derivatives is obtained. Because of the 
structure of the Kramers-Moyal expansion the block diagonal form of the matrix W in 
absence of the nonlinearity is expected to be unaffected by this deviation from white noise. 
Therefore we expect the main result of the work, namely the independence of the localization 
length of the nonlinearity to exhibit some degree of robustness and to hold for a wide range 
of models, beyond the specific model that is studied in detail in the present work. All these 
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problems should be subject of further studies. 
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